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ABSTRACT 

o ; 

Q , We determine and tabulate A^x\/Ak, the wavelength dependence of interstel- 

lar extinction, in the Galactic plane for 1.25/im < A < S.Oyum along two lines of 
CN ! sight: I = 42° and / = 284°. The first is a relatively quiescent and unremarkable 

pq ! region; the second contains the giant H II region RCW 49 as well as a "field" 

^ ! region unrelated to the cluster and nebulosity. Areas near these Galactic longi- 

1^ '. tudes were imaged at J, H, and K bands by 2MASS and at 3-8/im by Spitzer ioi 

the GLIMPSE Legacy program. We measure the mean values of the color excess 






O ! ratios {A[x] — Ak)/{Aj — A^) directly from the color distributions of observed 

Q ! stars. The extinction ratio between two of the filters, e.g. Aj/Ak, is required 

r-| ! to calculate Aix\/Ak from those measured ratios. We use the apparent JHK 

5^1 magnitudes of giant stars along our two sightlines, and fit the reddening as a 

O I function of magnitude (distance) to determine Aj/kpc, Ax/^pc, and Aj/Ak- 

^ [ Our values of A[x]/Ak show a flattening across the 3-8yum wavelength range, 

' roughly consistent with the Lutz et al. (1996) extinction measurements derived 
for the sightline toward the Galactic center. 

y^_ , Subject headings: dust, extinction - infrared: ISM 
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Introduction 



Extinction by interstellar dust affects most astronomical observations. The wavelength 
dependence of interstellar extinction has been studied extensively, but is not well-understood 
in the 3-9/im region (Draine 2003). The InfraRed Array Camera (IRAC; Fazio et al. 2004) 
on board the Spitzer Space Telescope is in the process of vastly increasing the number of 
observations in this wavelength region. Understanding the effects of dust extinction in the 
IRAC bands is important to properly interpret these observations. 

The wavelength dependence of interstellar extinction, Ax, is commonly treated as "uni- 
versal" in the infrared because it apparently varies far less between different sightlines than 
does extinction in the optical and ultraviolet. Many authors have concluded that Ax is a 
power-law {Ax oc A^'^) between ~l/im and ~4/im. Martin & Whittet (1990) found j3=1.8 in 
the diffuse interstellar medium (ISM) as well as the outer regions of the p Oph and Tr 14/16 
clouds. Other authors have fitted values between 1.6 and 1.8 (Draine 2003). A value of /3=1.8 
implies Ej_h / Eh_k={Aj — A^) / {Ah — Ak)=1-8±0.1 (the value changes slightly depending 
on the exact filters and source spectrum). Deep surveys of specific dark clouds in the JHK 
bands have revealed significantly different color ratios, ranging from Ej-h/Eh-k=^-4:7^0-0Q 
for luminous southern stars (He, Whittet, Kilkenny, & Spencer Jones 1995) to 2.08±0.03 
in the Coalsack (Racca, Gomez, & Kenyon 2002, and refs. therein). At wavelengths ap- 
proaching 9.7/im, Ax is dominated by absorption by the Si-0 stretching mode of interstellar 
silicates. There is uncertainty in the wavelength dependence of interstellar extinction be- 
tween the power-law regime at 1-2/im and the discrete silicate feature. Observations of H2 
rovibrational lines in Orion (Bertoldi et al. 1999; Rosenthal, Bertoldi, & Drapatz 2000) are 
consistent with continuation of the power-law to ^ 4/im, but their uncertainties are too large 
to constrain the extinction at longer wavelengths. Infrared Space Observatory (ISO) obser- 
vations of hydrogen recombination lines by Lutz et al. (1996) towards the Galactic center 
show a flattening of Ax in the region 3/im^ A ^dfira. Lutz (1999) confirms this fiattening 
using more recombination lines towards the Galactic center, and also shows evidence of extra 
mid-IR extinction towards external galaxies. 

The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; see Ben- 
jamin et al. 2003), a S'^iteer Legacy program, uses IRAC to obtain photometric images of the 
inner Galactic plane in four filters ([3.6], [4.5], [5.8], and [8.0]/im) simultaneously. By com- 
bining stellar photometry of the GLIMPSE images with the 2MASS^ point source catalog. 



^This publication makes use of data products from the Two Micron All Sky Survey (2MASS), which is a 
joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronautics and Space Administration and the National 
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we can sample A\/Ak at seven wavelength points [A] in the near-infrared (NIR, defined in 
this paper as 1.2-8.0;um): 1.240, 1.664, 2.164, 3.545, 4.442, 5.675, and 7.760/im. The adopted 
wavelengths are the isophotal wavelengths of the 2MASS and IRAC filters convolved with 
a K2III star, which we use as our standard probe of extinction along the line of sight (see 
discussion below) ^. 

We describe the data and their basic reduction in §2, and our method of measuring 
color excess ratios E[x]~k/Ej^k in §3. Section 4 discusses the conversion of these colors to 
Ax, and issues related to interpretation. In §5 we summarize and tabulate our recommended 
relative extinction values for the GLIMPSE/IRAC bands. 



2. Photometry of GLIMPSE data 

This paper uses stellar photometry in the Galactic plane obtained from two subsets of 
GLIMPSE data. The first dataset (Figure 1), 0.°33xl.°72 centered on (/,6) = (42°, 0.5°), was 
obtained during In-Orbit Checkout (IOC; October 2003, FID 631). The imaged region is 
a fairly unremarkable part of the Galactic plane. The second dataset (Figure 2), 283.9° < 
/ < 284.5° and -1.3° < b < 0.7°, was obtained for the purposes of Observing Strategy 
Validation (OSV; December 2003, FID 195). The region contains the giant H II region 
RCW 49 and its associated massive cluster (Churchwell et al. 2004, and references therein). 

Basic reduction of IRAC images was performed by the Spitzer Science Center (SSC) 
Pipeline (ver. 8.9.0 at / = 42° and ver. 9.0.1 at / = 284°). Positional accuracies are 
better than 1" (Werner et al. 2004). Point source full-width-half-max resolution of IRAC 
data ranges from ~1.6" at [3.6] to ~1.9" at [8.0]. The data were further processed by the 
GLIMPSE pipeline (Benjamin et al. 2003). Point sources were found and fiuxes extracted 
from each frame using a version of DAOPHOT (Stetson 1987), modified to be more robust 
in regions of rapidly varying high background (Babler et al., in preparation). To calculate 
stellar colors, we adopt zero magnitude fiux densities of 277.5, 179.5, 116.6, and 63.1 Jy for 
the four IRAC bands (M. Cohen 2004, private communication). For this present study, we 
culled the catalog to include only those sources with signal-to-noise greater than 10 in each 
band. The magnitude to which stars are well-measured in all bands is limited by the [5.8] 
and [8.0] bands that typically contain high diffuse backgrounds (the calculated uncertainty 



Science Foundation. Throughout this paper we wiU refer to the 2MASS Ks fiher as K. 

^Fiher wavelengths depend on the spectrum of the source, but differences are small and have little effect on 
our results - for example the flat-spectrum isophotal wavelengths for 2MASS are 1.235, 1.662, and 2.159/im, 
and the SSC-provided isophotal wavelengths for IRAC are 3.535, 4.502, 5.650, and 7.735/im. 
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Fig. 1. — GLIMPSE 3.6/im image of the region near / = 42° used in this study. The bright 
spots in the center are highly saturated stars; the regions around them are excluded from 
GLIMPSE photometry (marked circles). 
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Fig. 2. — GLIMPSE 3.6/im image of the region near / = 284° used in this study. The 
region near the giant H II region ROW 49 is marked with a circle, and the "field" region 
away from nebulosity, star formation, and significant CO emission (Dame, Hartmann, & 
Thaddeus 2001) is marked with a square. 



-6- 



for flux density is higher if the local background is higher). At the limiting magnitude of 
12.5, photometric uncertainties are less than 0.08, 0.09, 0.12, and 0.14 mag in the [3.6], [4.5], 
[5.8], and [8.0] bands, respectively. If we only require detections at wavelengths <5/im (e.g. 
for measuring the extinction at [4.5]), then we can include sources as faint as 14 mag at 
[3.6] and [4.5], with photometric uncertainties <0.15 mag. The flux calibration was checked 
by comparing extracted fluxes with modeled fluxes for flve early A-type dwarf stars - these 
agreed to within 7% for all IRAC bands (Cohen et al. 2003; Kurucz 1993). Detections of the 
same source in different images and different fllters were bandmerged (cross-identifled), and 
detections were additionally merged with the 2MASS point source catalog to provide JHK 
fluxes for many sources. Tests of bandmerging in crowded flelds using synthetic data show 
that fewer than 0.5% of sources are falsely cross-identifled down to the flux limits used in 
this work. 

The data segment at / = 42° was imaged before the telescope was fully focused. The 
point source function (PSF) is thus different from that in subsequent IRAC data. We do not 
expect the precision of our photometry to suffer since in the GLIMPSE pipeline the PSF 
is constructed from the data for each epoch of data collection. We did not assess the flux 
calibration independently for the IOC data as we did for the OSV data, but we do not see 
any gross errors, nor do we expect any, because the standard Spitzer calibration strategy 
was determined before launch. Results in this paper are derived from relative colors, which 
are insensitive to zero-point issues in the photometry. 



3. Measurement of Extinction 

We determine the relative extinction A^x]/^k using photometry of approximately 10,000 
stars ([A] refers to the three 2MASS and four IRAC fllters). Two different processes are re- 
quired. First, we measure color excess ratios E[x]-k/Ej-k by fltting the loci of a population 
of stars in color-color diagrams. Secondly, we need a single ratio of extinctions, say Aj/A^, 
determined or estimated independently. The calculation of A^x]/^k from E[x]-k/Ej-k fol- 
lows simply: 

Various observational studies and dust theories flnd a range of Aj/Ak between 2.25 and 
2.75 (Cardelli et al. 1989; Fitzpatrick 1999; Draine 2003). We determine the ratio Aj/Ak 
directly from our data by fltting the locus of red clump giants in a color-magnitude diagram. 
Both steps in our process are only possible using a fairly sensitive, uniform, large area survey 
such as 2MASS or GLIMPSE. 
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3.1. Determining Aj/Ak 

We determine the extinction ratio Aj/Ak by fitting the locus of reddened K giants in 
color-magnitude space. These stars are bright enough {Mj <^ —1.5, e.g. Lopez- Corredoira, 
Cabrera-Lavers, Garzon, & Hammersley 2002) that 2MASS and GLIMPSE can observe 
them to great distances through the Galactic plane. The intrinsic luminosity and color 
distributions of red clump K giants is determined by the evolution of stars from the main 
sequence and the initial mass function, is fairly narrow (± <^ 0.5 mag) and should not change 
with distance. We use the wavelength-independent effect of distance on apparent magnitude 
(the 1/d^ change in flux), together with the wavelength-dependent effect of extinction, to 
determine the average extinction per unit distance, and thus the absolute ratio of extinctions, 
Aj/Ak. 

In Figure 3, the apparent J magnitudes of the fleld stars near / = 284° are plotted 
against their J — K colors. This color-magnitude diagram shows features common along 
Galactic sightlines, in particular those at longitudes where the disk component dominates 
over the bulge component. A vertically extended concentration of points at J — K ^ 0.5 
are predominantly main-sequence dwarfs at a range of distances (8365 sources in this par- 
ticular plot). A second concentration (3560 sources) extending from J ~ 12, J — ii' ~ 1 to 
J~16,J — -ft'~2 are predominantly giants. Distance spreads the distribution vertically 
in the flgure; extinction moves a statistically uniform sample of stars vertically and hori- 
zontally by an amount proportional to distance, provided that the dust per unit distance is 
approximately constant. We note that the dwarf and giant sequences are not parallel - this 
is expected if the distribution of dwarfs is affected not only by interstellar reddening at faint 
J but also by the addition of intrinsically faint, numerous cool dwarfs. The absolute mag- 
nitudes of giants, by contrast, are dominated by the red clump giants of early K type, with 
fairly narrow color and luminosity distributions. There are no intrinsically faint local stars 
to be added to the giant sequence at faint J; instead, we have a population of intrinsically 
similar beacons, spread along the line of sight. The distinct red clump feature is discussed 
and exploited to probe Galactic structure by Lopez- Corredoira, Cabrera-Lavers, Garzon, 
& Hammersley (2002), and further used by Drimmel, Cabrera-Lavers, & Lopez-Corredoira 
(2003) to calculate the amount of extinction along certain Galactic lines of sight. The 859 
points redward of the giants in Figure 3 are probably mostly intrinsically red asymptotic 
giant branch stars. 

We perform a 4-variable minimization to flt the curved red giant branch star locus and 
determine the amount of extinction per unit distance. The curve can be written parametri- 
cally as a function of distance d: 

J = Jo + 51og(rf/10pc) + cj(rf/10pc) (2) 













8 


- 


■■■1 ■ 


■ ' '. , 


- 






" \' ■ '■ 


'* * * ■ 








' ' .\ *" 


^' 


- 






■ -''.V'" 


' ' ^« ' / ' ' " '. 


- 


1U 


- 




■t ■■,'.•■. ■ 


- 






•.■?;■• -^ 




- 








/^.■.■; •;.-■" .v '. ■ 


- 






t" 4 '"^^'t * * ' * 




12 


- 


m 


Kli.':'. . 


- 






w^ 


gS|?^^y >•>.. . . 


- 


14 


- 


s 


^^|¥: ■ . 


- 






^^^^^M 










'^H 










^ 


BhB^^-: 


- 


16 


(0) 









2 
J-K 



10 



-> 12 



14 



16 



(b) 







^ ?:?:-^v^- :;■■■■ 



2 
J-K 



Fig. 3. — The first panel (a) shows a color- magnitude diagram for field sources (away from 
RCW 49) near / = 284°. Features common in Galactic sightlines are the vertically extended 
concentrations of points at J-K ~0.5 corresponding to dwarf stars and at 1;^ J-K <^2 (red 
points) corresponding to mostly red clump giants (see text and Drimmel, Cabrera-Lavers, & 
Lopez- Corredoira 2003). We can isolate the red clump giants and fit their distribution (best 
fit is the overp lotted line) to determine Aj/A^ and the absolute extinction with distance (J 
magnitudes per kpc, see text for discussion of uncertainties). The second panel (b) shows 
the plot for the region around the giant H II region RCW49 - interestingly, the spread of 
red clump giants is very weak faintward of J ~14, J-K ^1.75. This is likely caused by an 
opaque interstellar cloud, and if that cloud is associated with RCW49, we can estimate the 
distance to be 5±1 kpc. 
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J-K = Jo-iro + (cj-c;,)(rf/10pc), (3) 

where cj is the average extinction per unit distance in the J band, assumed constant for a 
given hne of sight. We select the red giant clump stars (red points in Figure 3) by their JHK 
apparent magnitudes and colors, and use the IDL program amoeba to minimize the sum of 
the squared distances between each data point and the line defined by Equation 3. If we 
wish only to determine the ratio cj/ck = Aj/Ak, then Jq and Kq can be determined by best 
fit for the population of giants along the particular sightline. We get the same results for 
cj/ck whether Jq and Kq are fixed at the values suggested by Lopez-Corredoira, Cabrera- 
Lavers, Garzon, & Hammersley (2002) for a K2III {M^ = —1.65, H ~ K = 0.75, intrinsic 
spread of 0.3 magnitudes in the population), or allowed to vary. We varied the region of the 
sky, the color selection criteria for RGB stars, the initial guess for the four parameters, and 
the degree of robustness in the fit algorithm (iteratively rejecting high sigma outliers). The 
range in fit coefficients is larger than the numerically calculated sigma for each fit, and we 
use that larger range as our quoted uncertainty. The uncertainty in the extinction per unit 
distance cj is larger than in the ratio cj/ck, because it depends on assuming the absolute J 
and K magnitudes of the red clump. 

We fit the locus of red clump giants in all possible JHK color-magnitude diagrams (J — 
K versus J, H—K versus K, etc). Our derived ratios are Aj/Ak='2.5±0.2, Ah/Ak=1.55±0A, 
and Aj/Ah=1.Q5±0A. As noted above, we also measure the average magnitudes of extinc- 
tion per kpc along this line of sight: cj=0.35±0.15, c/f=0.25±0.1, and c/^=0.15±0.1 magkpc"^. 
These are about la higher than numbers often used for the average Galactic disk, cy ~0.7 mag kpc" 
or cj ~0.2 mag kpc~^. Examining the color-magnitude diagrams also reveals two regions 
of high extinction along the line of sight towards / = 284°: at J ~12.7, J — K ~1.3 and 
J ~14.3, J — K ~1.7. The extinction of the more distant cloud is so high towards the 
ROW 49 giant H II region that few red clump giants are seen beyond it (see Figure 3). If 
that cloud is associated with ROW 49, as seems reasonable, then we can estimate the dis- 
tance to be 5±1 kpc, independently of any other distance measurement. Finally, it is more 
difficult to fit the giant branch towards / = 42° because of greater contamination by AGB 
stars and the bulge red giants at magnitudes fainter than J ~ 14, but we find the same 
ratios of Aj/Ah/Ak within uncertainties. 



3.2. Determining color excess ratios 

We measure each color excess ratio E[x]~k/Ej^k by fitting the locus of a population of 
stars in a color-color diagram. The color excess E[x]^k is defined as the difference between 
a star's observed and intrinsic color. The excess for a particular star is proportional to 
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the amount of dust along the sightline to that star, but ratios of color excesses reveal the 
wavelength dependence of extinction. If one plots colors [A] — K versus J — K, the slope of 
the source distribution is the color excess ratio E[x\-k/Ej-k, as illustrated in Figure 4. 

In each region of the sky, we fit a line to the stellar colors ([A] — K)/{J — K) of only 
the giant stars (selected by their JHK magnitudes and plotted in red in Fig. 3). Extremely 
red sources (at the bottom of the first 4 panels in Figure 4) likely have nonphotospheric 
emission from circumstellar dust; these outliers must be excluded from measurement of 
interstellar dust extinction, because the dust emission has a different wavelength dependence 
and appears as nonphysical extra reddening. A small number of highly evolved (AGB) 
stars can also fall in the excluded region because of absorption bands in their atmospheres. 
We calculated fits using several different algorithms, and used the differences in results to 
estimate our uncertainties; the simplest fit performed was a linear regression weighted by 
uncertainties in both axes. In addition to fitting all sources and fitting only red clump giants, 
we tried arbitrarily excluding sources with extreme red excesses (e.g. with ([A] — K) < 
— 0.3(J — K) — 0.65). We also used a robust algorithm which iterates the fit, rejecting 
high-sigma outliers. This also has the effect of excluding some of the extreme red excess 
sources. 

Figure 5 shows the measured color excess ratios E[x]~k/Ej-k plotted as a function of 
A~^ for the RCW 49 giant H II region, the field near / = 284°, and the field near / = 42°. 
The plotted values are listed in Table 1, along with average values. The uncertainties are 
a combination of variation in fit results depending on the exact part of the sky, stellar 
population, and fit algorithm, calculated uncertainty of the fit coefficients, and errors in 
GLIMPSE and 2MASS photometry (in order of decreasing significance). There is good 
agreement among the excess ratios measured in all three regions, although the data are from 
very different directions in the Galactic plane. More extensive interpretation is possible after 
conversion of the measured color excess ratios into relative extinction A[x]/^k, as discussed 
in the next section. 

An alternative method for measuring the color excess ratios due to interstellar extinction 
that is not dependent on the intrinsic stellar color distribution is described by Kenyon, Lada, 
& Barsony (1998). The method consists of comparing a selection of stars presumed to be 
mostly behind an extincting cloud with a selection of stars not extinguished by the cloud 
("off-cloud"). Instead of plotting colors (e.g., H — K), one plots all possible color differences 
([A] — K)i — ([A] — K)j versus {H — K)i — {H — K)j (where i indexes all "cloud" stars and 
j all "off-cloud" stars). If the stellar color distributions are statistically the same in the 
two samples, the color difference plot is a tight distribution extended along the extinction 
vector. As before, the slope of that distribution is the color excess ratio E[x\^k/Ej_k- We 
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Fig. 4.— Color-color diagrams of sources detected by GLIMPSE and 2MASS in the "field" (6 > 0, 
away from known H II regions) near I = 284°. The slope of a line fitted to each source distribution 
is the color excess ratio E[x]-k/Ej~k- We fit only the red clump giants, selected by their location 
in the J, J-K color-magnitude diagram and plotted here in red. The first four panels (a)-(d) show 
source colors J-K versus [X\-K for the four IRAC bands [A]. Panels (d) and (e) show J-K versus 
H-K, for sources with high-quality fiuxes in the JHK bands, and for only those sources with 
high-quality fluxes in all 7 filters, respectively. If one fits only the latter stars, the slopes change by 
^5% and this has been included in error estimates. The dominant source of scatter in these plots is 
variation in the intrinsic colors of stars (~ ±0.3 mag Lopez-Corredoira, Cabrera-Lavers, Garzon, & 
Hammersley 2002). Sources with nonphotospheric emission from circumstellar dust (traditionally 
called "infrared excess" sources) appear as extreme outliers, e.g. at the bottom of the first panel 
\H.O]-K versus J - K. 



-12- 



calculated the wavelength dependence of extinction from GLIMPSE and 2MASS data using 
the RCW 49 region as "cloud" and the region away from the H II region to the north as 
"off-cloud." The resuhs {E[x]^k/Ej^k = 0.36±0.02, -0.31±0.02, -0.36±0.02, -0.40±0.03, 
and -0.37±0.02, for H, [3.6], [4.5], [5.8], and [8.0], respectively) are the same within the 
line-fitting uncertainties to the values listed in Table 1 that were obtained from the simpler 
method described above. 

We have discussed ratios of color excesses E[x\~k/Ej^k- The use of J and K in the de- 
nominator is convenient because we determine Aj/A^ from the color- magnitude relation (see 
§3.1), but the A[x\/Ak we derive is independent of this choice. For completeness, we fit other 
color excess ratios, and found that the resulting ratios did not change within uncertainty - 
for example 



El4.5]-K 
Ej-K 
-^[4.5]- [3.6] , -,\ f E[3,Q]^K 



0.35 ±0.02, 



E[3.6]-K J V Ej^K 

(0.16 ± 0.03 + 1)(0.30± 0.02) = 0.35 ± 0.07 



4. The wavelength dependence of interstellar extinction Ax 

4.1. Determination of Ayx^/Ax 

With color excess ratios from the slopes of color-color plots (§3.2) and Aj/Ak deter- 
mined from the color-magnitude relation of giants (§3.1), Equation 1 yields A\^x^/ Ak- The 
presence of enough data points from 2MASS and GLIMPSE to measure Aj/A^ directly 
from our data is very important to this study. Color excess ratios are sensitive only to 
changes in the extinction with wavelength. Without Aj/Ak-, we would have to extrapolate 
E[x]-k/Ej-k in Figure 5 from K band to A~^ = 0; the intercept of that extrapolated curve 
is —Ak/{Aj — Ak) = 1/(1 — Aj/Ak)- Although this extrapolation is a commonly used 
procedure, it is difficult to perform accurately, and questionable because of the absorption 
features of silicates peaking at 9.7 and 18/im. The arrowhead at the lower left of Figure 5 
shows the value of 1/(1 — Aj/Ak) using our values measured as described in section 3.1. By 
assuming a power-law through the JHK relative extinctions (effectively an extrapolation 
to A^^ = 0), Martin & Whittet (1990) estimated Aj/Ak = 2.7. This value requires the 
intercept in Figure 5 to be -0.59. Our intercept derived from fitting the red clump giants is 
at -0.67. 

Figure 6 shows the derived A[x]/Ak for the field stars at / = 284° (solid diamonds). 
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Fig. 5. — The color excess ratio E\^x\_k/ Ej^k plotted as a function of 1/A, for the RCW 49 
region (dot-dashed) and for the "field" near / = 284° (solid), and for the "field" near / = 42° 
(dashed). Error bars on the "field" curves are smaller than for RCW 49, and comparable to 
the sizes of the points. Horizontal bars indicate the filter widths and the fact that we have 
only sampled the wavelength dependence of interstellar extinction at those discrete filters. 
The filled arrowhead is discussed later in the text. 



-14- 



for / = 42° (open diamonds), and for the dust surrounding RCW 49 (solid circles). The 
uncertainties are a combination of uncertainty in determining Aj/A^, variation in fit results 
depending on the exact part of the sky, stellar population, and fit algorithm, calculated 
uncertainty of the fit coefficients, and errors in GLIMPSE and 2MASS photometry (in order 
of decreasing significance). There is good agreement between the three determinations of 
A[x]/Ak, in two very different directions of the Galactic plane, and even towards a giant H II 
region. 

A simple fit to the average A[x]/Ak in the three regions is given by 

log [Ax/Ak] = 0.61(±0.04) - 2.22(±0.17) log(A) + 1.21(±0.23)(log(A))^ (4) 

where A is in /xm. The fit curve is plotted along with the data points in Figure 6 and is 
recommended for use between 1.25 and 7.75/xm. (Note that a simple power-law extinction 
curve with /3=1.8 has coefficients of 0.60, -1.8, and 0.0). 

For historical reasons, the filter to which Ax is compared is usually taken to be Ay, and 
our extinction measurements are presented as A^xi/^k- In contrast to NIR extinction, optical 
extinction (A < 1/im) is known to vary significantly among sightlines (Cardelli et al. 1989; 
Fitzpatrick 1999); the variation is well-characterized by a single parameter, Ry = Ay/Es^v- 
The variation of Ay /A^ is beyond the scope of this paper, but Ay/A^ ~8.8 for i?y=3.1 and 
~7.5 for Ry='b (Cardelli et al. 1989). In general, the lack of variation of the extinction curve 
throughout the NIR among various environments (in comparison to its strong dependence 
on Ry in the optical/ultraviolet) suggests that it may not be possible to estimate Ay /A^ 
for our dataset from NIR observations alone. 



4.2. Comparison with previous measurements and discussion 

Previous measurements of the wavelength and spatial dependence of interstellar extinc- 
tion in this wavelength regime differ. It is important to note that the measured excesses at 
IRAC wavelengths (A > 3/im) lie above the curve extrapolated from A < 3;um(see Figure 5). 
This implies a source of extinction in the IRAC wavelength regime that is in excess of the 
power-law wavelength dependence exhibited by a grain model tuned to UV-optical extinction 
observed in the diffuse ISM (e.g. Weingartner & Draine 2001, Case A). Our measurements 
reject a pure power-law extinction curve by >4o" at 6/im. Lutz et al. (1996) and Lutz (1999) 
observed the Galactic Center with the Short Wavelength Spectrometer on the Infrared Space 
Observatory (ISO) in the interval 2.4-45/im, using recombination lines of hydrogen to pro- 
vide intrinsic flux ratios. They converted their reddening values to Ax/Ak by assuming a 
theoretical value of A^/ Ay from Draine (1989). Figure 6 shows that the Lutz et al. (1996) 
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Fig. 6. — In the first panel, symbols show the derived relative extinction values for RCW 49 
(solid circles), the field near / = 284° (solid diamonds), and the field near / = 42° (open 
diamonds). The open and solid diamonds coincide at A=3.5/im. The lines show fits (see 
Equation 4) to the mean field (solid) and RCW 49 (dot-dashed). Horizontal bars indicate 
the filter widths and the fact that we have only sampled the wavelength dependence of 
interstellar extinction at those discrete filters. In the second panel, the measurements of Lutz 
et al. (1996) have been added (squares with error bars), as well as theoretical curves from 
Weingartner & Draine (2001) for i?v'=5.5 (their case "A" dashed, their case "B" dotted). 
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Table 1. IRAC Extinction A[a]M 



ii^ 





J 


H 


K 


3.6 


4.5; 


5.8; 


8.0 




1.240/im^ 


1.664^m 


2.164/im 


3.545/im 


4.442/im 


5.675/im 


7.760/im 




E[x]-k/Ej. 


-K 












RCW49 


1.0 


0.37±0.02 


0.0 


0.31±0.03 


0.38±.03 


0.40±0.05 


0.37±0.04 


/ = 284° 


1.0 


0.35±0.02 


0.0 


0.30±0.02 


0.35±.02 


0.38±0.02 


0.38±0.03 


/ = 42° 


1.0 


0.36±0.02 


0.0 


0.29±0.02 


0.37±.02 


0.35±0.03 


0.35±0.04 


average 


1.0 


0.36±0.02 


0.0 


0.30±0.01 


0.37±.02 


0.38±0.04 


0.37±0.04 




A[x]/Ak 


RCW49 


2.50±0.15 


1.56±0.06 


1.00 


0.53±0.07 


0.42±0.09 


0.40±0.14 


0.45±0.11 


I = 284° 


2.50±0.15 


1.54±0.05 


1.00 


0.57±0.05 


0.43±0.07 


0.41±0.07 


0.37±0.07 


/ = 42° 


2.50±0.15 


1.54±0.06 


1.00 


0.57±0.06 


0.45±0.07 


0.48±0.09 


0.48±0.11 


average 


2.50±0.15 


1.55±0.08 


1.00 


0.56±0.06 


0.43±0.08 


0.43±0.10 


0.43±0.10 



'^The adopted wavelengths are the isophotal wavelengths of the 2MASS and IRAC filters 
convolved with a K2III star - see 51. 
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values are very similar to those presented in this work but show more scatter. 

Our values of A[x]/Ax are consistent with Lutz et al. (1996) and Lutz (1999) in spite of 
the large differences in method and sightlines. Their color excess ratios pertain to the line 
of sight to the Galactic center, while our / = 42° sightline probes only the outer 30% of that 
galactocentric distance and / = 284° is confined to approximately the solar circle. The path 
to the Galactic center contains Ay ~ 5 — 10 magnitudes of extinction from molecular clouds 
(which may exhibit different extinction than the diffuse ISM; see discussion and references 
in Whittet 1997), out of a total Ay ~ 25 magnitudes. The molecular fraction of the ISM 
in our regions is at present unknown at adequate spatial resolution; low- resolution ^^CO 
(Dame, Hartmann, & Thaddeus 2001) measurements suggest that only a few magnitudes of 
the observed Ay ~ 10 — 15 may be from molecular clouds. Finally, using nebular emission 
lines to provide the extinction assumes that the extinction is uniform over the 14" x 20" 
aperture of the ISO spectrometer, whereas using stellar photometry probes a pencil beam 
along each line of sight. Our measurements are consistent with Lutz et al. (1996) despite all 
of these differences in methodology and circumstance. 

The similarities of the wavelength dependence of extinction for our two very different 
Galactic longitudes and, to within larger errors, towards the Galactic center (Lutz et al. 1996; 
Lutz 1999), is remarkable. This suggests an almost universal extinction law in the infrared. 
Even towards the massive star-forming region ROW 49 the extinction is very similar to that 
observed in the "field" regions. 

It is useful to estimate the extent to which the [8.0] IRAC band is affected by the 
silicate feature that peaks at 9.7/im. Most stars emit the peak of their emission far short- 
ward of the IRAC wavelengths, and their emission in the Rayleigh- Jeans limit weights the 
filter response towards shorter effective wavelengths. The Spitzer Science Center-provided 
isophotal wavelength is 7.735/im, at which the silicate opacity is <10% of that at 9.7/im (e.g. 
Jaeger et al. 1994, 1998). The ISO spectra of heavily obscured objects in Gibb, Whittet, 
Boogert, & Tielens (2004) show 8/im fluxes that seem little affected by silicate absorption. 
Nevertheless, the filter is broad, and silicate absorption in the diffuse ISM appears consistent 
with a FWHM of 2.3yum, possibly wider in molecular clouds (Roche & Aitken 1984; Bowey, 
Adamson, & Whittet 1998). Convolution of the IRAC filter profile with a source emitting a 
Rayleigh- Jeans spectrum and with a Gaussian absorption feature at 9.7/im shows that the 
silicate feature could have as much as a 20% effect on the filter flux. Studies of individual 
objects (especially those with circumstellar dust) could show more variation in the [8.0] ex- 
tinction with source spectral shape than the average ISM results presented here. The effect 
of the silicate feature on the [5.8] IRAC band is negligible. 

The ratio of polarization to extinction is an important diagnostic of the nature of grains. 
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Consider aligned grains of a particular type of material. At a particular wavelength, both 
extinction and polarization are integrals of similar optical constants over the size distribution 
of the grains. Both polarizations and extinctions at both 4 and 8/im have been measured in 
two molecular clouds: GL2591 and the BN object in Orion (see Martin & Whittet 1990 for 
observations and references.) Because of atmospheric opacity, there are no measurements 
between these wavelengths. In GL 2591, p{4 fim)/p{8fira) ~2.5; in BN, p{4fim) /p{8fim) ~2.0. 
The silicate strengths are also similar; p{10fim)/p{8fim) is ~3.7 for GL2591 and 4.3 for BN. 
Equation (4) shows that our measured extinction ratio A{4fim) / A{8fira) ~ 1.32, significantly 
below the 4yum/8/im polarization ratio of either object. The wavelength dependence of the 
polarization outside of the ice band is similar for these two clouds but differs significantly from 
our measured wavelength dependence of extinction, suggesting that the grains providing the 
8 iim extinction are different from the silicates at lOjim and those providing ^iim extinction 
(carbonaceous?), and that they are not as well aligned. The alternative to this conclusion 
is that the dust within Orion BN and GL 2591 has a wavelength dependence of extinction 
significantly different from the dust observed by GLIMPSE. 



5. Conclusions 

We have measured the wavelength dependence of interstellar extinction in the 1.25- 
8.0/im region by combining GLIMPSE and 2MASS observations. The values calculated for 
two very different directions in the Galactic plane are remarkably similar. Even near the 
giant H II region ROW 49, extinction is not significantly different. It may be possible to 
study variations of extinction with Galactic longitude using the entire GLIMPSE survey, 
expected to be complete by early in the year 2005. Since the values of Ej_h/Eh-k in the 
Galactic plane are somewhat higher than at higher Galactic latitude, there may be differences 
in the extinction at IRAC wavelengths as well. 

The similar derived behavior of the extinction along several different sightlines suggests 
that our relative extinction values A^xi/^k may be generally applicable. Our average extinc- 
tion measurements may be used to correct GLIMPSE and other IRAC data in the Galactic 
plane. These values are provided in Table 1, and the simple formula of Equation 4 provides 
an analytical fit of simple functional form. 
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